MPI-PhT/97-21 
FSU-SCRI-97-24 



Further results on O(a) improved lattice QCD 
to one-loop order of perturbation theory 

Stefan Sint 

SCRI, The Florida State University, 
Tallahassee 32306-4052, Florida, U.S.A. 

Peter Weisz 

Max-Planck-Institut fur Physik 
Fohringer Ring 6, D-80805 Miinchen, Germany 

Abstract 

We present results at one-loop order of perturbation theory for various 
improvement coefficients in on-shell 0(a) improved lattice QCD. In particular 
we determine the additive counterterm required for on-shell improvement of 
the isovector vector current. Employing a general mass-independent renor- 
malization scheme we also obtain the coefficients of the 0(a) counterterms 
which are proportional to the quark mass in the improved isovector pseudo- 
scalar, axial vector and vector operators. In the latter case a comparison with 
a recent non-perturbative study is made. 
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1. Introduction 



This paper belongs to a series of publications on chiral symmetry and 
0(a) improvement in lattice QCD with Wilson quarks [1-6]. Here we report 
on the computation of various improvement coefficients to one-loop order of 
perturbation theory. The underlying theoretical framework is based on the 
Schrodinger functional (SF) [7,8] and has been presented in detail in the papers 
referred to above. In particular, a general overview can be obtained from 
refs. [1,2]. 

The present paper is organized as follows: In section 2 we recall various 
definitions limiting ourselves to those which are essential for the understanding 
of this paper. These include in particular the definitions of the on-shell 0(a) 
improved isovector pseudo-scalar, axial vector and vector operators, and of 
the renormalized parameters and fields in a mass-independent renormalization 
scheme. We then summarize our main results which are contained in table 1 
and eq. (2.16). 

In section 3 we present a few details of the one- loop calculation and discuss 
the renormalization procedure for the general case of non-vanishing renormal- 
ized quark mass. Section 4 contains the determination of the improved vector 
current and of various improvement coefficients which arise in the case of non- 
zero quark mass. We end with a few concluding remarks (section 5). Finally, 
two appendices have been included which provide analytic expressions for the 
tree-level amplitudes and counterterms at zero and non-zero quark mass re- 
spectively. 



2. Overview and results 

A general introduction to on-shell 0(a) improvement of lattice QCD with 
Wilson quarks can be found in ref. [3]. In particular there it has been discussed 
how the renormalization procedure can be carried out in a way consistent with 
0(a) improvement. In the following we adopt the conventions and notations 
as introduced in this reference. 
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2.1 Definitions 

On-shell 0(a) improvement requires the introduction of 0(a) counterterms 
for both the lattice action and the composite fields of interest. The improved 
lattice action can be obtained by including the Sheikholeslami-Wohlert inter- 
action term (with coefficient c sw ) [9]. Concerning the composite fields we will 
restrict attention to a few gauge invariant combinations which are bilinear in 
the quark fields. Assuming Nf > 2 mass degenerate quark flavours we consider 
the local isovector fields t, 

V^x)=i;(xh^T a ^x), (2.1) 

A*(x) =^(x) 7M 7 5 |r a V(x), (2.2) 

P a {x) =^(x) 75 ir a V(x), (2.3) 

T^(x) =i^(x)a^r a iP(x). (2.4) 

Using these definitions and following ref. [3] the improved isovector vector and 
axial vector currents may be parametrized as follows, 

(Vi)l =V; + c v (gl)a\{dl + d v )T£ v , (2.5) 

(Aj)l =A% + c A {gl)a\{dl + d,)P a . (2.6) 

Throughout the paper we will use a mass-independent renormalization 
scheme which is compatible with 0(a) improvement. In such a scheme the 
renormalized coupling and quark mass are given by [3] 

9i=giZ s {gl,afi), (2.7) 

m R =rh q Z m (gQ,afi), (2.8) 

where \x is the renormalization scale. The parameters go and m q are defined 
as follows 

g%=go[l + b e (g%)am q \, (2.9) 
rh q =m q [l + M#o) am q] > m q = m o - m c (2.10) 

t The Dirac gamma matrix conventions are as in appendix A of ref. [3] 
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where m c is the critical bare quark mass. 

With these definitions the renormalized improved currents and pseudo- 
scalar density take the form [3], 

(V R )l =Z v (glafi) [1 + by{gl)am^] (Vft, (2.11) 

(A R )° =Z A (glafi)[l + b A (g 2 )am (l ](A l )«, (2.12) 

(P R ) a =Z P (g 2 ,a f i)[l + b P (g 2 )am cl ]P a . (2.13) 

In the massless theory on-shell correlation functions involving only these 
fields are 0(a) improved if the improvement coefficients c sw , ca and cy are 
properly chosen as functions of the bare coupling g . In the presence of massive 
quarks one also needs to know the various 6-coefficients introduced above. 

To compute these coefficients we will consider correlation functions de- 
rived from the Schrddinger functional which involve the above operators and 
also the boundary quark fields (,( and (',(' [7]. 0(a) improved correlation 
functions are then obtained with the renormalized fields Cr.;Cr an d Cr>Cr> 
which are all related to the bare fields by the same renormalization factor, e.g. 

Cr = Z c (g 2 , a/i) [1 + b c (g 2 )am q ] (. (2.14) 

For completeness we recall that Schrddinger functional boundary condi- 
tions lead to additional cutoff effects which can be cancelled by counterterms 
that are localized at the boundaries. These boundary counterterms can be 
chosen such that the associated improvement coefficients c Sjt and c s ,t appear 
as weight factors for certain terms of the lattice action close to the bound- 
aries [7,3]. 

2.2 Results 

The computation of c sw to one-loop order of perturbation theory was first per- 
formed in ref. [10] and the numerical result was confirmed to three significant 
digits in ref. [4]. For the purpose of our one- loop calculation it is sufficient to 
know that c sw = 1 at tree-level of perturbation theory. 
The improvement coefficients ca and cy are given by 

ca = c^g 2 + O(^), 4 1} = -0.005680(2) x C F , (2.15) 
cv = c^g 2 + O(ocJ), 4 1} = -0.01225(1) x C F , (2.16) 
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Table 1. Tree level and one- loop improvement coefficients b 



X 


b (0) 

u x 


°X 


g 





0.012000(2) x Nf 


m 


1 

2 


-0.07217(2) x C F 


c 


1 
2 


-0.06738(4) x C F 


A 


1 


0.11414(4) x C F 


V 


1 


0.11492(4) x C F 


P 


1 


0.11484(2) x C F 


T 


1 





where C F = (A 2 — 1)/2N for gauge group SU(A). The coefficient has 
first been obtained in ref. [4] and the computation of Cy- 1 will be presented 
in section 4. The extended numerical data produced in the course of this 
calculation also led to the improved estimate for as given in eq. (2.15). 
For the various ^-coefficients, whose perturbation expansion reads 

b = b (0) +b W g 2 + O{g ^ (2 . 17) 

we have collected the known results in table 1. The computation of b^ was 
performed in ref. [11] and the calculation of the remaining coefficients is the 
topic of section 4. It is remarkable that the coefficients bty and bp '' are so 
close numerically. To our knowledge there is, however, no obvious theoretical 
explanation of this fact. Note also that all coefficients appear to be reasonably 
small, i.e. of 0(1) when expanded in powers of <?o/47r for gauge group SU(3). 

In the case of by we can now make a comparison with the non-perturbative 
result for N{ = and gauge group SU(3) which was obtained in ref. [6]. An 
updated fit of the data, which reduces to the correct one-loop expression for 
small go, is given by T 

h - 1 ~ °- 7613g ° 2 + - 0012g ° 4 ~ °- 1136g ° 6 < a < 1 (2 18) 

bv( 9o )- 1 _ o.9145 5o 2 ' °^o<l. (2.18) 



We thank Hartmut Wittig for providing this fit and fig. 1. 
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Fig. 1. The improvement coefficient by as a function of the bare 
coupling for gauge group SU(3). The data points and the fit function 
corresponding to the dashed curve are taken from ref. [6]. The solid curve 
is the new fit function (2.18), the dotted line is the perturbative result and 
the crosses denote the perturbative result when the coupling gp is used 
instead of go. 

The data together with the fit functions are displayed in fig. 1. Note that both 
curves are indistinguishable for all practical purposes. 

Using the one-loop result for by we may also compare the data with 
so-called boosted perturbation theory. To this order this simply amounts to 
replacing the bare coupling g by Parisi's coupling g p [12], defined by 

gl = g 2 /P, (2.19) 

where P denotes the expectation value of the plaquette in infinite volume. The 
corresponding values of 6y are plotted as crosses in fig. 1. Although the use 
of g P moves the perturbative result towards the non-perturbative curve, there 
remains a significant discrepancy at the larger values of the bare coupling. 
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3. The one-loop calculation 



We define the relevant correlation functions and provide some details of 
their evaluation to one-loop order of perturbation theory. This is followed 
by a discussion of the renormalization procedure for the general case of non- 
vanishing renormalized quark mass. In the following the reader is assumed to 
be familiar with ref. [4]. In particular the notation and the general set-up of 
perturbation theory on a finite lattice with Schrodinger functional boundary 
conditions will be taken over from this reference. 

3.1 Definition of the correlation functions 

We consider the following correlation functions which were introduced in 
refs. [1,3], 

/a(* ) = -a 6 ^i(^(x)C(y)7 5 ir a C(z)), (3.1) 
fp{x ) = -a 6 ^i(P a (x)C(y) 7 5ir a C(z)). (3.2) 

Here and in the following we adopt the convention that repeated indices are 
summed over. To study the improved vector current (2.5) we also introduce 
the new correlation functions, 

k Y (x ) = -a 6 ^i(V^(x)C(y) 7 4^CH), (3.3) 
y,z 

k T (x ) = -a 6 ^i(T fc a (x)C(y)74^C(^))- (3-4) 
y,z 

The amplitudes above are sufficient for the determination of most of the im- 
provement coefficients of sect. 2. However, in order to carry out a few addi- 
tional checks we also calculated the boundary-to-boundary correlation fi , 

h = ~lJ £ ^(C / (u)7 5 ir a C , (v)C(y)75ir a C(z)), (3.5) 

u,v,y,z 

which has previously appeared in the normalization conditions for the isovector 
axial vector and vector currents in refs. [1,6]. As an aside we recall that /i 
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will also be needed for the computation of the running quark mass in the 
SF scheme as defined in ref. [1]. A detailed discussion of this topic will be 
presented elsewhere [13,14]. 

3. 2 Integration over the quark fields 

The fermionic action being bilinear in the quark fields, the corresponding 
Grassmann integration can be carried out analytically using Wick's theorem. 
To write down the resulting expressions in a compact form we follow ref. [4] 
and introduce the matrix H(x) through 

Here ip c i denotes the classical solution of the Dirac equation in a given gauge 
field configuration and p is its boundary value at time xq = 0. One then finds 

fA(x ) = -±(tr{H(x)^ H(x)}) G , (3.7) 
f P (x ) = l(tT{H(x)^H(x)}) G1 (3.8) 

kv(x ) = -i<tr{ 75 7fc^(x) t 757fc^(^}> G , (3.9) 
k T (x ) = i(tr{ 757fc #(x)t 757fc7o tf(x)}} G , (3.10) 

where the trace is over the Dirac and colour indices. The bracket (• • -)g means 
that expectation values have to be taken with the effective gauge field measure 
including the fermionic determinant. 

In order to obtain a compact expression for fi we follow ref. [6] and 
introduce the matrix 

Li — ) xo=T — a 

x 

Then one has 

h = \{tr{K^K}) G . (3.12) 
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3.3 Perturbation expansion 

The perturbation expansion is now easily generated following ref. [4]. In par- 
ticular we choose vanishing boundary gauge fields C and C and take over the 
corresponding gauge fixing procedure. The gluon field is introduced in the 
standard way by parameterizing the link variables according to 

U(x,fi) = exp{agoqfj,(x)}. (3.13) 

As explained in ref. [4], the classical quark field and thus H(x) can be expanded 
in perturbation theory, 

H(x) = H^(x)+g H^(x)+gf ) H^(x) + O(g 3 ). (3.14) 

Also expanding the boundary improvement coefficient [3], 

c t = l + c^g 2 + O(gt), (3.15) 

the corresponding expansion for K reads 

K = KW+goKW+fiKW + O( 5o 3 ), (3.16) 

with 

^ (0) = $E p +^ (0) ( T - a > x )> ( 3 - 17 ) 

X 

3 

RW = 73 Y, p +\ Hil Hx) ~ ago(x)tf (0) (*)) , (3.18) 

_L — I ) xo=T — a 

x 

R{2) = p E p + { R{2) (*) - a *> ^ R{1) ( x ) 

X 

+ ^[aq (x)] 2 H^(x) + c { t 1) H^(x)} . (3.19) 

> xo=T — a 

Inserting these expansions in eqs. (3.7)-(3.10),(3.12) and carrying out the in- 
tegrations over the gluon and ghost field variables finally leads to the desired 
expansion for the correlation functions, 

/ = / (0) +5o 2 / (1) + 0(5 4 ), (3.20) 
where / stands for any of the amplitudes in eqs. (3.1)-(3.5). 
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3.4 Tree level results 

Explicit expressions for the tree level correlation functions and f p ^ have 
been given in ref. [4]. The study of their approach to the continuum limit led 
to the determination of b^,b m ,bx, bp and ca at lowest order of perturbation 
theory. To this order we have the relations 

4%o) = |/ P 0) (*o) - (3.21) 

4 0) (xo) = f/ A 0) (xo)-i^ 0) K), (3.22) 

A (0) = {|/r^o)-i/f(xo)} , (3.23) 

> x Q =T — a 

and we then find the tree level results by^ = b^ = 1 and Cy' 1 = 0. All 
coefficients are thus known to lowest order of perturbation theory (cf. sect. 2). 

3.5 Computation of k^ and 

Except for the different Dirac structure the calculation of ky^ is completely 
analogous to the cases of and fp^ treated in ref. [4]. In particular, there 
are again three diagrams to be computed which are the same as in fig. 1 of 
this reference (where the cross in this case denotes the vector vertex). The 
contribution of the quark boundary counterterm can also be inferred from 
there by noting 

k { y ] (x ) h = |/ P 1) (x )b - |/i 1) (x )b. (3.24) 

Note that the correlation function k^ could be computed along the same 
lines. However here, &t is only encountered as an 0(a) counterterm for the 
correlation function of the improved vector current (2.5). Since the associated 
improvement coefficient cy vanishes at tree level we will only need the lowest 
order expression k^ in the following. 

The correlation function fi differs from the cases treated previously and 
therefore deserves a more detailed presentation. Inserting the expansion of 
the matrix K (3.16) in eq. (3.12) we obtain 

fV = \{tr{K^KV}) & 

+ \{ tv{K^K^ + K^K^})^ (3.25) 
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x =T 



x =T-a 




x =0 11 1 1 

3 4 5a 5b 



x =T 
x =T-a 



x =0 




6a 6b 7a 7b 



Fig. 2. Diagrams contributing to f± at order g%. The dotted lines 
symbolize the link from Euclidean time T — a to T. 

Here, the bracket {■ ■ - )q means the integration over the gluon and ghost fields 
as explained in ref. [4]. While the integration over the ghost fields is trivial at 
this order of perturbation theory, the integration over the gluon fields generates 
the diagrams displayed in fig. 2. The first three diagrams are similar to the 
ones computed for /a,/p and ky, and the explicit link variables in eq. (3.11) 
lead to the additional diagrams 4-7. 

Furthermore, receives a contribution from the quark boundary coun- 
terterm which can be written in the form (see also eq. (4.31) of ref. [4]) 

fS = 2c1 1} A ( ° } + hfp \ T ~ «)b " ~ a) b . (3.26) 
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An explicit evaluation shows that this expression is indeed of order a, as 
expected. 

On a lattice of a given size T/a x (L/a) 3 the Feynman diagrams can be 
evaluated numerically by inserting the explicit time-momentum representation 
of the propagators and vertices into the expressions for each diagram. For this 
purpose two independent Fortran programs were written and the final results 
were checked against each other and for gauge invariance. A fast version 
of one of the Fortran programs then enabled us to obtain numerical results 
for a relatively large range of lattice sizes, thus allowing for rather precise 
extrapolations to the continuum limit (cf. sect. 4). 

3.6 Renormalized correlation functions 

In order to take the limit of large L/a at fixed ratio T/L and physical length 
L, the bare parameters and fields have to be scaled such that the renormalized 
parameters and fields stay fixed. Furthermore this has to be done in a way 
consistent with 0(a) improvement. It has been shown in ref. [3] that these 
requirements are met by any mass-independent renormalization scheme with 
the properties as summarized in sect. 2. In particular, the 0(a) improvement 
coefficients introduced there are then independent of the renormalization con- 
ditions. 

In perturbation theory a particularly convenient choice is the minimal 
subtraction scheme on the lattice which we shall adopt in the following. In 
this scheme all renormalization constants at a given order in the perturbation 
expansion in eft are polynomials in ln(a/x) with mass-independent coefficients 
and no constant parts. 

To first order of perturbation theory the substitutions for the coupling 
constant and the quark mass then amount to 



9o=9i + 0(gi), 



(3.27) 



mo = m^ +9 imP + 0( 9 i), 



(3.28) 



where the precise form of the coefficients 




(3.29) 



m. 



(i) 



(i) _ 




(3.30) 



o 



i/l — 2am R 



11 



is a direct consequence of the definitions made in sect. 2. 
The renormalized correlation functions, defined by 



[M*o)]r = Z Y (1 + 6 v aw q )Z c 2 (l + b^arriq) 



x {ky(x ) + ac v \(d*Q + d )k T (xo)} , 



(3.31) 



[fp(xo)]n = Z P (l + & P am q )Z 2 (l + o c am q ) 2 jp(x ) 
[/a(^o)]r = ^a(1 + b A am q )Z% (1 + 6 c am q ) 2 



(3.32) 



x {/a(zo) + ac A |(<9o + d )fp(x )} , 



(3.33) 



[/ 1 ] R = Z c 4 (l + 6 c am q ) 4 /i, 



(3.34) 



have a well-defined perturbation expansion in the renormalized coupling <7 R , 
with coefficients that are computable functions of a/L. However, to determine 
the 0(a) improvement coefficients we only need to know the renormalized 
amplitudes up to terms of order a 2 . Neglecting such terms, the expansion of 
[ky] R reads 



where it is understood that all quantities on the r.h.s. are evaluated at mo = 
mj ' as given in eq. (3.29). Analogous expressions are obtained for [/a] r and 
[/p]r> i n particular the case m R = has already been discussed in ref. [4]. We 
thus directly proceed to the result for [/i] R , 



[ky(x )] R 





(3.35) 



[ /i]r = (1 - 2am R ) [/<°> + + m^£- o f[ 0) 

+ (4Z« + am R [4bf + 2Z«]) A (0) }] + 0(g R ). (3.36) 
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Since we are neglecting terms of order a 2 , the expansion 



rn. 




tor + am n 



(z«+6£>)]+0(a a ) 



(3.37) 



may be inserted in eqs. (3. 35), (3. 36). However, note that we do not directly 
expand the coefficient TOq " 1 [eq. (3.29)]. One has to be careful here because the 
bare one- loop amplitudes are linearly divergent. While the additive renormal- 
ization of the quark mass removes these divergences, the correct evaluation of 
the renormalized amplitudes to 0(a) requires a consistent treatment of TOq - 1 
to order a 2 , in both the bare one-loop amplitudes and the mass counterterms. 
It thus appears safer to first carry out the renormalization procedure using 
the exact coefficient and only neglect terms of 0(a 2 ) in the final result for the 
renormalized correlation functions. 



The improvement coefficients can be determined by requiring an improved 
continuum limit behaviour of the renormalized correlation functions. We first 
consider the case of vanishing renormalized quark mass and compute Cy ^ . We 
then turn to the massive case and also determine the 6-coefficients. 

4-1 Choice of parameters 

For the numerical analysis of the renormalized amplitudes one needs to make 
a choice for the kinematical parameters T,L,xq,6 and the quark mass m R . 
Following ref. [4] we choose T = 2L and xq = T/2. Setting z = m R L we define 
the dimensionless functions, 



4. Computation of the improvement coefficients 



h P (9,z,a/L) = [f P (x )] R 



x =T/2> 



(4.1) 



h A (6,z,a/L) = [/aOeo)] 



R 



x =T/2> 



(4.2) 



h v (9,z,a/L) = [k v (x )] R 



x =T/2 



(4.3) 



h dA (6,z,a/L) = L±(d% + d )[f A (x )] 



R 



xo=T/2> 



(4.4) 



fei(e,z,a/L) = [Mr. 



(4.5) 
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From the discussion in sect. 3 one then infers the form of hp, h^A and hi 
hp = u + 5r|""i + c[^U2 + arnipuz 

+ (z£> + 2Zf + + 2b™]) u } + 0(4), (4.6) 

^dA = Wq + 4 j w i + c{^w 2 + am^w 3 + c^w^ 

+ (Z« + 2Zf + ±z[b<» + 26«]) «,„} + 0(4), (4.7) 

hi = t + 4 j*i + c| 1) t 2 + am^h 3 

+ (4zf) + |.[4&« + 2Z«]) to} + 0(4). (4.8) 
Similarly we find, 

= «o + 4 j^i + Ct 1} f2 + am^v 3 + c^V 

+ (z£> + 2Z« + + 2bf}) vo} + 0(4), (4.9) 

and the same for hy, with the subscript A replaced by V and the coefficients 
Vi replaced by yj. 

All the coefficients are still functions of 6 and z. Analytic expressions can 
be derived for those coefficients which derive from the tree level correlation 
functions or the 0(a) counterterms. Their expansions to order a/L for the 
special case of z = and non-vanishing 6 are given in appendix A and in 
appendix B of ref. [4]. In appendix B we give the corresponding results for 
8 = and non- vanishing renormalized quark mass. 

The coefficients vi,yi,ui,wi and t\ are only obtained numerically so that 
we must choose specific values for the parameters 6 and z. In the case of 
vanishing quark mass z = 0we decided to collect numerical data for the three 
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values 6 = 0, 6 = 0.1 and 6 = 1. In the massive case we set 6 = and chose 
the three values z = 0.1, 2 = 0.5 and z = 1. 

With the above choices for the parameters the Feynman diagrams were 
then evaluated numerically in 128 bit precision arithmetic for a sequence of 
lattice sizes ranging from Lj a = 4 to L / a = 48. 

4-2 Logarithmic divergences 

In the minimal subtraction scheme the one-loop renormalization constants all 
take the form 

Z(glafi) = 1 + Z^f Q + O^ 4 ), (4.10) 

where Z^> is proportional to \n(afi). From analytical results obtained with 
conventional perturbative techniques one then expects [15-17,8] 

4 1} = 4 1} = o, (4.ii) 

and 

Z^ = -Z^ = -2Z^ = ^Ha,). (4.12) 

In the course of our computation we were able to determine the coefficients 
of the logarithmic divergences numerically to high precision. The results were 
in complete agreement with the above expectations and will not be further 
discussed. In the numerical analysis presented below we set fi = 1/L and use 
the exact coefficients when subtracting the logarithmically divergent parts. 

4-3 Computation of cty 

For vanishing quark mass the analogue of eq. (4.9) for hy simplifies to 
hy(6,0,a/L) = y + g\ 



2/i + c ( ^y2 + am { pyz 



+ (4 1 ) + 2Zf ) ) m + 4 1} y 4 ] + 0{gi). (4.13) 



In order to extract Cy^ one can for example demand O(a) improvement of the 
combination hy — |(2/ip — /ia)- Due to the identities (3.21) and (3.24) the 
coefficients and omf do not contribute and the improvement condition 
then reads 

?/i - |(2m - Vi) - lZ^u - Key - c { l ] )v A = const + 0(a 2 ). (4.14) 
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Here we have anticipated the absence of 0(a) effects proportional to ln(L/a), 
which are cancelled by setting the coefficient of the Sheikholeslami-Wohlert 
term to c sw = 1. We have checked in a number of cases that this indeed 
happens. Defining 

C(L) = yi - i(2«i - Vl ) + ^u ln(L/a), (4.15) 

and the symmetric difference operator 

\{d + d*)C(L) = i- \C{L + a) - C(L - a)] , (4.16) 

we find 

<#> -<£>=- lim ^- {d + d * )C (L). (4.17) 

It turns out that the 0(a/L) corrections in the expression above are large 
for both values of 9 considered and the extrapolation must therefore be done 
carefully. We used the method described in ref. [18], appropriately adapted to 
the present case and finally obtained 

C W - c ^ 1} = -0.00657(1) x C F . (4.18) 

ic 

quoted in eq. (2.16). 

An alternative determination of Cy" 1 can be obtained by requiring 0(a) 
improvement of the ratio hy(0, 0, a/L)/hy(0, 0, a/L). With the results of ap- 
pendix A this is equivalent to demanding 

Vi- ^yo(yi) e =o + am i 1] (ys + 2yo^) +c ( v ) y 4 + c { t 1) y 2 = const + 0(a 2 ). (4.19) 

While the renormalization constants cancel in this ratio we now need to know 
the one-loop coefficient mf of the critical bare quark mass. We here follow 
ref. [4] and use the relation 

(w 1 /w 3 ) e=z=0 = -am,™ + 0(a 3 /L 3 ). (4.20) 

We then define the combination 

C'(L) = Vl - ly ( yi ) e=0 - (y 3 + 2y ^) ( Wl /w 3 ) e=z=0 , (4.21) 



Insertion of the numerical value for [eq. (2.15)] then leads to the result 
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where it is understood that w\ and u>3 are computed at the same lattice size 
as the other coefficients. With the results of appendix A we then find 



" 4 C ° sh(2 ^>: 1 af) = li,n ^( fl + ,T,m> 

COsh(2V3 6') i/^oo 1/4 



Evaluation of this relation at two different values of 9 allows to eliminate 
and thus leads to a direct determination of Cy ^ . Complete consistency with 
the previous method was found. 

With our extended data we now also obtain a slightly more accurate 
estimate of than the one given in ref. [4], viz. 

c£ 1} = -0.01346(1) x C F . (4.23) 



4-4 Computation of the b- coefficients to one-loop order 

In order to compute the 6-coefficients to one-loop order we now consider the 
massive case where z = m-^L is kept fixed at a non-zero value. In the following 
we always assume 9 = and use the coefficients given in appendix B. 

We start by considering hi [eq. (4.8)], which can be written in the form 



h x = Ne 



-4z 



1 + 



gl{T(L) + 4z(b™ + ct } + + 0(4)] • (4-24) 



Using the exact results for the logarithmic divergences we have defined 



T(L) = N- 1 ^ 



ti + am^H^ 



+ 2(l + 2*)^ln(L/a). (4.25) 



To evaluate T{L) numerically on a lattice of size L/a, we here again follow the 
procedure above and replace the coefficient am?' by the ratio — {wi/w^)e= z =o^ 
computed for the same lattice size. Requiring improvement of hi then deter- 
mines the linear combination of improvement coefficients, 



4z 



c (1) 



+ bf ) +zb^ 



lim (L 2 /a)Ud + d*)T(L). 



L/a — > 00 



(4.26) 



Since c[^ is already known, b^ and b^ can be determined by evaluating the 
improvement condition at two different values of the quark mass. 
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We can proceed similarly for the other correlation functions fp , /a and 
ky. We here require the improved continuum limit behaviour for the ratios 
hp /hi, h\/hi, and h\jh\. Defining 



U(L) = N^e 22 



V(L) = -N^e 22 
y(L) = N~ l e 2z 



vi + am^V3 



+ 2z J^HL/a), 



(4.27) 



+ (l + 2z)^ln(L/a), (4.28) 



+ (l + 2z)^ln(L/a), (4.29) 



we find e.g. 



1 h 



v 



1 + gl 2y(L)-T(L) + 2zy(24 ) +b^) +0(4). (4.30) 



N hi 



L 



Similar expressions can be derived for the other cases and we thus find that 
0(a) improvement requires 

zb p 1] = lim {L 2 /a)\{d + d*) \U{L) - \T(L)} , (4.31) 

L/a— >oo L J 



2c£ ) +# ) = lim (L 2 /a)±(d + d*) }>(L) - ±T(L) . (4.33) 



lim (L 2 /a)I(a + 9*) V(L) - \T{L) , (4.32) 

L/a— >oo L J 



L/a— >oo 



The evaluation of the various improvement conditions was done using the 
methods of ref. [18] and led to the numerical results as quoted in table 1. 
Various consistency checks were made by studying different combinations of 
correlation functions. For example, a more direct determination of was 
obtained by considering the ratio hdA(0,z,a/L)/h,A(0,z,a/L). As a further 
confirmation of the conceptual framework, we have also analyzed data at finite 
renormalized quark mass and values of / and verified that the ^-coefficients 
are indeed independent of 0. 
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5. Concluding remarks 



This paper presents a further step in a systematic investigation of the 
continuum limit of lattice QCD. At one-loop order of perturbation theory we 
have studied the continuum limit of many different renormalized amplitudes. 
Various non-trivial consistency checks have been carried out, for both zero and 
non-zero renormalized quark mass. As a result our analysis further supports 
the consistency of the general framework of on-shell 0(a) improvement as 
presented in ref. [3]. 

The 0(a) improvement coefficients which enter the definition of the im- 
proved lattice action and the quark bilinear isovector pseudo-scalar, vector 
and axial vector operators are now all known to one-loop order of perturba- 
tion theory. The ultimate aim is of course to compute the same coefficients 
non-perturbatively. This has been partially achieved for the quenched theory 
(i.e. for Nf = 0) in refs. [5,6]. More precisely, the coefficients c sw ,ca and 6y 
have been determined for this case and a similar determination of cy is under 
way [19]. The extension to the full theory with two quark flavours has also 
been initiated and first results are expected in the near future [20]. 

However, the non-perturbative methods developed so far cannot be used 
in a straightforward manner to determine all of the improvement coefficients. 
An example is the coefficient 6a for which one presently has to rely on per- 
turbative estimates. Of course, this does not represent a severe limitation 
provided the quark mass (measured in lattice units) is small enough. In this 
context it may be of some interest that the perturbative results for 6y and 6a 
are approximately of the same magnitude. While this need not remain true 
at higher orders or beyond perturbation theory it is certainly tempting to use 
the non-perturbative result for 6y as a first guess for 6a- 
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Appendix A 



For the massless case z = m K L = 0, the expansions of the functions Ui,Vi 
and Wi (for i 7^ 1) have been given in appendix B of ref. [4], up to corrections 
of order (a/L) 2 . Our definitions then imply 

Vi = - \vi for i = 0, 2, 3, (A.l) 
Vi = -\vi + §u>o-^- (A.2) 
Using the abbreviations 

co = cosh(2v / 3 6 | ), si = sinh(2\/3 6), (A.3) 
we also give the explicit expressions, 



iY ( 2CO+1 ) (KA\ 

Vo = \ 2 , A.4 
3co^ 



8iV0si(co + l) a , k x 

^ = I' (A - 5) 



iVsi(co + 2)L r 7V0si(co + 2) 
3a/35co 3 a I lSy^co 3 



+ — ^(co 3 + 4co 2 - 2co-6H-, (A.6) 
2N6si a , A . 

* " TwZ' (A ' 7) 



For the functions ij introduced in eq. (4.8) we obtain 



to = (A.8) 



^ _ 24iV6»si a 
\/3co 3 

2iVsi L 2 AT 
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(A.9) 



cl9N0si 8N8 2 (l-2si 2 )ia 

13^3 W (A - 10) 



Appendix B 

In this appendix we provide explicit expressions for the coefficient func- 
tions introduced in eqs. (4.6)-(4.9), for the special case 6 = 0. We set 
z = m K L and use the tree level expression of the renormalized quark mass, 
m K = mo(l — amo/2). Up to terms of order a 2 /L 2 we then find 



t = Ne~ iz 



u 3 



t 3 = Ne~ 4z {-4^ + 2 + Az + 2z{\z 2 + 3z - l)|}, 



Uo = Ne~ 2z 



u 2 = Ne~ 2z 2zY, 



Ne- 2z {-2^ + 2z + z 2 (lz-l)^-}, 



-Ui for i = 0, 2, 3, 



t, 4 = _JVe 2z 2zy, 
u, = iVe -22 2z, 



W3 = Ne~ 2z { (2 - Az) ^ + 4z 2 - 2z 
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B.l 
B.2 

B.3 

B.4 
B.5 

B.6 

B.7 
B.8 

B.9 



w 2 = Ne~ 2z Az 2 y, (B.10 



+ . 2 (5-fz+|. 2 )|} 



(B.ll) 



W4 = Ne 




(B.12) 



for 



i = 0,2,3,4. 



(B.13) 
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